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Abstract 

The search for a canonical set of eigenvectors of the discrete Fourier transform has been ongoing 
for more than three decades. The goal is to find an orthogonal basis of eigenvectors which would 
approximate Hermite functions - the eigenfunctions of the continuous Fourier transform. This 
eigenbasis should also have some degree of analytical tractability and should allow for efficient 
numerical computations. In this paper we provide a partial solution to these problems. First, we 
construct an explicit basis of (non-orthogonal) eigenvectors of the discrete Fourier transform, thus 
extending the results of [7]. Applying the Gramm-Schmidt orthogonalization procedure we obtain 
an orthogonal eigenbasis of the discrete Fourier transform. We prove that the first eight eigenvectors 
converge to the corresponding Hermite functions, and we conjecture that this convergence result 
remains true for all eigenvectors. 
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1 Introduction 

The fractional Fourier transform is becoming increasingly more important due to an ever-growing list of 
applications in signal processing, optics and quantum mechanics and in other areas of Science [1, 13]. In 
order to define this object one needs first to diagonalize the Fourier transform operator V, and that is 
where Hermite functions 


d n 

bnP) := (— n n\)- 1/2 e x2/2 —e- x \ n > 0 (1) 

dx n 

become indispensable. It is well-known that Hermite functions are the eigenfunctions of the Fourier 
transform operator 

{Ffi n ){x) \= —E= f e~ lxy fi n (y)6.y = (-i) n fi n (x), 
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dirn(V) 
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L 

L 
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L + l 

L 

N = AL + 3 

L + l 

L + l 

L + l 

L 


Table 1: Dimensions of the eigenspaces {W, X, Y, Z} of the discrete Fourier transform Xn, corresponding 
to the eigenvalues {1, —i, —1, i}. 


and that they form a complete orthonormal basis of L 2 (M). In other words, the orthogonal basis of 
Hermite functions diagonalizes the Fourier transform operator X and this allows us to define the fractional 
power (-F) 0, - the fractional Fourier transform. Of course, there are infinitely many ways to choose an 
eigenbasis of the Fourier transform, each of them would lead to a different version of the fractional 
Fourier transform. What makes Hermite functions so special among all other eigenfunctions is that 
they are analytically tractable and they have many useful properties. This analytical tractability and 
Mehler’s formula [1] (which gives an explicit expression for the integral kernel of the fractional Fourier 
transform) have sealed the standing of Hermite functions as the canonical set of the eigenfunctions for 
Fourier transform. 

Compared to the continuous Fourier transform, the situation with the eigenvectors of the discrete 
Fourier transform (DFT) is much more complex. Let us first introduce several necessary definitions and 
notations. The (centered) discrete Fourier transform Tn is a linear map that sends a vector u G into 
a vector v e according to the rule 

«(() = -j= Y. e- 2 ’ m/N u(k), l 6 I N , (2) 

kei N 


where In denotes the set of integers 

In ■= {—M + 1, —M + 2,..., —M + N}, M= [(N + l)/2\. 

Everywhere in this paper we will follow the convention that all vectors will be denoted by bold lower 
case letters and the elements of a vector will be labeled by the set In- Before we begin discussing the 
eigenvectors of the DFT, let us say a few words about its eigenvalues. The eigenvalues of both continuous 
and discrete Fourier transform are given by ±1, ±i. In the continuous case each eigenspace is infinite 
dimensional, whereas the multiplicities of eigenvalues in the discrete case are presented in the Table 1. 
This result was obtained in 1921 by Schur, however it can be easily derived from a much earlier result of 
Gauss on the law of quadratic reciprocity, which is essentially a statement about the trace of the matrix 
corresponding to the linear transformation Xn■ The formula for the multiplicities of eigenvalues was 
rediscovered by McClellan and Parks [9] in 1972, and since then there have appeared many other proofs 
(see [12] for a very simple proof based on the Vandermonde determinant formula). 

Returning to our discussion of eigenvectors of the DFT, it seems that the first example of an eigenbasis 
was constructed by McClellan and Parks [9] as a by-product of their proof of the results in Table 
1. Other constructions of eigenvectors of the DFT include a number-theoretic construction [11] and a 
representation theory approach [6, 15]. Mehta [10] has constructed an eigenbasis based on theta functions. 
One important approach for finding orthogonal eigenvectors of the DFT is based on certain symmetric 
tridiagonal (or almost tridiagonal) matrices which commute with the DFT. The eigenvectors of these 
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matrices give an orthogonal eigenbasis of the DFT. Such matrices were first discovered by Dickinson and 
Steiglitz [4] and by Grunbaum [5] (see also [2, 3, 14] for more recent developments). 

The above discussion clearly shows that there is an abundance of different sets of eigenvectors of 
the DFT. How can we choose “the best” among them? First we need to define what we mean by “the 
best”. An ideal set of eigenvectors of the discrete Fourier transform J-'n should satisfy the following 
three requirements: the eigenvectors should be orthogonal, they should converge to the corresponding 
Hermite functions as N —> +oo and they should have at least some degree of analytical tractability 
(so that we are able to compute them). Judging by these standards, most of the above candidates 
can be disqualified. For example, the eigenbasis of McClellan and Parks is explicit but not orthogonal 
and it does not approximate Hermite functions; the eigenbasis of Mehta [10] does approximate Hermite 
functions but it is not orthogonal; the eigenbasis constructed in [15] is orthogonal, but it requires N to 
be a certain prime number and it does not seem to approximate Hermite functions. In fact, the only 
sets of eigenvectors that do satisfy the above requirements come from the commuting symmetric matrix 
approach of Dickinson and Steiglitz [4] and Grunbaum [5] , though in both cases the eigenvectors are not 
given explicitly and to find them one needs to diagonalize a tridiagonal symmetric matrix. 

Not much is known about explicit Hermite-type eigenvectors of the DFT. In fact, it seems that the 
only result in this direction comes from the recent paper by Kong [7], who gives two examples of such 
eigenvectors. More precisely, when N = AL + 1 the vector u G with 

2 L 

u{k) = (cos(27r/c/iV) — cos(27tj /N )), k G In, (3) 

j=L +1 

is an eigenvector of Fn, and the same is true when N = 4L for the vector v G with 

2L—i 

v(k) = sin(27r k/N) (cos(27r/c/iV) — cos(27t j/N)) , k G In- (4) 

j=L+ 1 

Numerical evidence given in [7] suggests that these two vectors approximate the first two Hermite func¬ 
tions. What is also interesting about these two eigenvectors is that they are unique in a certain sense. 
It follows from the proof of the main result in [7] (though it is not explicitly stated there) that the only 
eigenvector of the DFT of the dimension N = 4L+1 which has zero elements u(k) = 0 for L+l < \k\ < 2L 
is the vector u defined in (3). 

In this paper we provide the first known example of an explicit orthogonal eigenbasis of the DFT 
for which the individual eigenvectors seem to approximate the corresponding Hermite functions. First 
of all, we construct a basis of which consists of even and odd vectors similar to Kong’s vectors u 
and v given above. These vectors behave nicely under the DFT and we prove that they are unique in 
a certain sense. The proof of these results uses the methods of Kong [7] combined with the g-Binomial 
Theorem. From this basis of we construct an eigenbasis of the DFT following the same procedure as 
McClellan and Parks [9]. As a corollary of this construction, we obtain a new proof of the result on the 
multiplicities of the eigenvalues of the DFT presented in Table 1. Finally, we apply the Gramm-Schmidt 
orthogonalization procedure to the basis of each eigenspace of the DFT. This gives us an orthonormal 
basis of , consisting of the eigenvectors of the discrete Fourier transform Fn- We prove that the first 
eight eigenvectors converge to the corresponding Hermite functions as N —y +oo and we conjecture that 
this convergence holds true for all eigenvectors. 

The paper is organized as follows. In Section 2 we present our results for the case N = l(mod 4). 
The proofs of these results are collected in Section 3 and in the Appendix we give the corresponding 
formulas for all remaining cases N G {0,2,3}(mod 4). 


3 


2 Results 


We begin by stating some definitions and notations that will be used throughout this paper. We denote by 
W, X, Y, Z C R N the eigenspaces of the DFT corresponding to eigenvalues 1, —i, —1, i (their dimensions 
are given in Table 1). The dot product of two vectors u and v will be denoted by u • v and the Euclidean 
norm of a vector in C N is defined as ||w || 2 — u ■ u. We denote by [xj and \x\ the floor and the ceiling 
function. From now on we will drop the subscript N in the notation of the discrete Fourier transform and 
will write simply T instead of Tn - there should be no confusion as the continuous Fourier transform 
will not be used anymore. 

We say that a vector a G R N is even (odd) if a(k) = a(j ) (respectively a(k ) = — a(j )) for all k,j G In 
such that k = — j (mod N). An equivalent definition can be given as follows: when N = 2 M + 1 the 
even (odd) vectors satisfy the condition a(k) = a(—k) (respectively, a(k) = —a(k )) for —M < k < M. 
When N = 2 M the even vectors satisfy a(k) = a(—k) for — M + 1 < k < M — 1 (note that there is 
no restriction on u(M )) and the the odd vectors satisfy a(k) = —a(—k ) for — M + 1 < k < M — 1 and 
a(M ) = 0. The following two properties are well known and will be used extensively later: if u G R N is 
an even vector and v G R N is an odd vector, then 

(i) Ju (i.Pv) is a real even (respectively, real odd) vector; 

(i) W 2 u = u and F 2 v = —v; 

The following object will appear frequently in this paper: the sequence {S(k)}k>o is defined by 

k 

S(k) f° r ^ (5) 

3 =1 

and S'(O) := 0. This sequence satisfies 

(i) S(k) = 0 for k > N; ( 6 ) 

(ii) S(k)S(N - k - 1) = N for 0 < k < N - 1. (7) 

Property ( 6 ) follows at once from the definition (5) while property (7) can be derived by writing the 
polynomial P(w) = (1 — w N )/( 1 — w) as a product of — 1 linear factors and then taking the limit as 
w —> 1 . 

We denote by supp(b) the support of a vector b, which is the set of all indices k G In such that 
b(k ) 7 ^ 0. Following Kong [7] we define the signal length (or simply, the length ) of a vector b as 

Z(b) := max{7 — j : i,j G supp(b)}. 

We record the following important relationship between the length of a vector and its DFT: For any 
nonzero vectors b G CW we have 

l(b)+l(Pb) > N + l. (8) 

This result (which reminds one of the uncertainty principle) follows from Theorem 1 in [7]. We include 
the proof of this result in Section 3 for convenience of the reader. 

The next theorem is our first main result. It gives the exact number of even/odd vectors that are 
“extreme” in the sense that the quantity l(h) + l(J-"b) is the smallest possible among all even/odd vectors. 
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Theorem 1. Assume that N > 3. 

(%) If N = 2M + 1 there exist only M + 1 nonzero even vectors satisfying /(u) + /(-Tu) — N + 1 and 
only M nonzero odd vectors such that /(v) + / (-Tv) < N + 3. 

(ii) If N = 2 M there exist only M — 1 nonzero even vectors satisfying u(M ) = (-Tu)(M) = 0 and 
/(u) + /(.Tu) < N + 2 and only M — 1 nonzero odd vectors such that /(v) + / (-Tv) < iV + 2. 

The results and definitions that we have presented so far are valid for all integer values of N > 1. 
From now on we will concentrate on the case when N = 1 (mod 4). The corresponding results for all 
remaining cases N G {0, 2, 3} (mod 4) are collected in the Appendix. 

2.1 The case when N = AL 1 

Our next goal is to find explicit expression for the vectors described in Theorem 1 and to compute their 
discrete Fourier transform. 

Theorem 2. 

(i) The even vectors {u n }-L<n<L described in Theorem l(i) are given by 

u n (k) S(3L + n + k)S(3L + n — k) for —L<n<L and k G In- (9) 

These vectors are linearly independent and they satisfy 

-Tu n = N~ 1/2 S(2L + 2n)u_ n . (10) 

(ii) The odd vectors { v n }-L<n<L-i described in Theorem l(i) are given by 

v n (k) := sin(27rfc/A^)S'(3L + n + k)S(3L + n — k) for — L < n < L — 1 and k G In- (11) 
These vectors are linearly independent and they satisfy 

-Tv n = -±N- l/2 S{2L + 2 n + l)v_ n _ 1 . (12) 


Remark 1 . Using the sum-to-product formula for the cosine function and identity (7) it is easy to check 
that formula (9) is equivalent to 


2 L 

u n {k) = N2 L+n (cos(2nk/N) - cos(2nj/N)). (13) 

j=L-n -\-1 


In particular, the vector u 0 is a scalar multiple of the vector u discovered by Kong [7] (see (3)). 

Let us highlight some important properties of the vectors u n and v„. First of all, property (6) implies 
that 


supp(u„) = {k G Z : \k\ < L — n} and supp(v T1 .) = supp(u n ) \ {0}. (14) 
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(a) 


(b) 


Figure 1: The red (respectively green, blue) circles correspond to (a) u n (fc)/||u n || (b) u n (A:)/||v n || 
for n — 0 (respectively n = —5, n — 5) and the black curve is the graph of the function (a) 
x (2/iV) 1 / 4 exp(— ttx 2 /N) and (b) x (->• 2 5 / 4 7T 1 / 2 lV~ 3 / 4 xexp(— nx 2 /N). Here IV = 41 and the x-axis 
represents k G In- 


In particular, we have Z(u n ) = /(v re ) = 2 L — 2n + 1 and it is clear that the vectors u ra and v n satisfy 

l(lln) + u n) — N + 1, I(V„) + l(J- V n ) = N + 3, 

so these are precisely the vectors described in Theorem 1. The vectors u n and v n should be viewed as 
the discrete counterparts of the Gaussian functions exp(—ax 2 ) and xexp(—ax 2 ) (see Figure 4). There 
are indeed many analogies. First of all, the elements of these vectors, Uo(k) and Vo(k), converge to the 
corresponding Gausian functions as N —y +cx) (this follows from Lemma 2 in Section 3). Second, the 
identities ( 10 ) and ( 12 ) are discrete analogues of the integral formulas 

[ e~ ax2 x e~ lzx dx = s/^fa x e -" 2/(4a) , f xe~ ax2 x e^dx = -1^/(4 a 3 ) x ze~ z2/{4a) . 

J ]R J M 

Finally, the vectors u n and v n are connected by the following identity 

u n (k + 1) — u n {k — 1) = — 4 sin( 7 r( 2 L + 2n) /N)v n -i(k), for — L < n < L and —2 L < k < 2 L, (15) 

which we will establish in Section 3. This result is the finite difference counterpart of the formula 

d~ ax2 = -2 axe~ ax2 . 
dx 

Our next goal is to construct a basis for each eigenspace W, X, Y and Z of the DFT. This construction 
is based on Theorem 2 and the following property: if u is an even vector (v is an odd vector) then J-'uiu 
(respectively, L7-V ± v) is an eigenvector of the discrete Fourier transform T with the corresponding 
eigenvalue ±1 (respectively, =Fi). This property can be easily established using the fact that T 2 u = u 
and X 2 v = — v for every even vector u and every odd vector v. 
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(a) n = 0 (black) and n = 1 (red) (b) n = 2 (black) and n = 3 (red) 



-20 -15 -10 -5 0 5 10 15 20 


(c) n = 4 (black) and n = 5 (red) 


Figure 2: Comparing the first six eigenvectors <f n (k) (black circles and red squares) with the scaled 
Hermite functions e l ^ 2 if n {ek) (black and red lines). Here N = 41, e = \J2jtJN and the x-axis represents 
the index k G [-20,20]. 


Corollary 1. Define N vectors {w n } 0 <n<L, {x n }o<n<L-i, {yn}o<n<i-i, {zn}o<n<L-i as foliouis: 

w n = u n + N~ 1 ^ 2 S(2L + 2n)u_ n , 0 < n < L, 

x„ = v n + N~ 1/2 S(2L + 2 n + l)v_ n _i , 0 < n < L — 1, 

y n = -u n+ i + N~ 1/2 S(2L + 2n + 2)u_ n _! , 0 < n < L - 1, (16) 

z n = — v n + iV _1 / 2 S'(2L + 2n + l)v_ ri _ 1 , 0 < n < L — 1. 


The vectors {w n }o<n<L (respectively, {x n } 0 < n <L_i ? {y„}o<n<L-i and {z n } 0 < n <L_i y ) form a basis for the 
eigenspace W (respectively, X, Y and Z). 

Remark 2. The results of Corollary 1 imply that the eigenvalues {1, —i, —1, i} of the discrete Fourier 
transform have corresponding multiplicities {L + 1 ,L,L,L}, which agrees with Table 1. 

The eigenbasis constructed in Proposition 1 is not orthogonal. Our next step is to build an orthonor¬ 
mal basis via the Gramm-Schmidt algorithm. 

Definition 1. Define N vectors {0fc}o<fc<iv-i as follows: 

_ w n _ X n ^ _ y n ^ _ Z n (-\>7\ 

Tin ~ |i ~ |i ) W4n+1 “ || ~ || ) W4n+2 — n ~ n ) W4n+3 ~ \\~ n ) GO 

11 ^n || || || 11 y n \ \ 11 11 

where the vectors {w n } 0 <n<L, {xJo<n<L-i, {y n }o< n <L-i, {zn}o<n<L-i are obtained by Gramm-Schmidt 
algorithm starting from {w„} 0 <n<L, {x„} 0 <n<L-i, {y n }o<n<L-i, {z n }o<n<L-i: 


w m = w m 


x„ = x n 


m —1 



w j, 0 < n < 


0 < n < L — 


L, 

1 , 
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Z„ = z n - " J Zj, 0 < n < L - 1. 

i=o H Zj 'H 


The following theorem is our main result in this paper. Recall that Hermite functions ip n (x) are 
defined by (1). 

Theorem 3. 

(i) The vectors {0 n } o < n <jv-i form an orthonormal basis in R N consisting of eigenvectors of the discrete 
Fourier transform: for all 0 < n < N — 1 we have = (—i ) n <f> n 

(ii) For 0 < n < 7 we have 


max 

—2L<k<2L 


e 1/2 (j) n (k) - if n (ek)\ ->• 0 as N —> Too, 


(18) 


where e := a JFkJN . 

The first six (respectively, twelve) eigenvectors <p n and the corresponding scaled Hermite functions 
are shown on Figure 2 for N = 41 (respectively, Figure 3 for N = 401). The results presented on Figure 
3 suggest that the convergence in (18) is not restricted to the first eight eigenvectors, and that it should 
hold for all of them. Thus we formulate 

Conjecture 1: For all n > 0 we have 

max I e~ 1 ^ 2 cb n (k) — tb n (ek)\ —>• 0 as IV —» Too. (19) 

—2L<k<2L 

Our proof of Theorem 3 - which is presented in Section 3 - is essentially a proof by verification, and 
it does not provide an intuitive reason as to why this convergence should hold. Our approach would 
probably work for n > 8, though larger values of n would require a prohibitive amount of algebraic 
computations. Thus our method of proof of Theorem 3 is probably not the right way for proving 
Conjecture 1. We hope that some other properties of the vectors u n and v n (which we may have 
overlooked) will lead to a simple and insightful proof of this conjecture. 

Figure 4 provides a “bird’s eye view” of the eigenbasis <f n and compares it with Hermite functions. 
In particular, on Figure 4c we show the values of n and k where 4> n (k) is either positive, or negative or 
equal to zero. From these results it seems that the vectors (p n have the same number of zero crossings 
as the corresponding Hermite functions ip n . Let us make this statement precise. We define the number 
of zero crossings of a vector a e 1 N as the number of pairs (i,j) such that (i) i,j G In, (ii) i < j, 

(iii) a(i)a(j ) < 0 and (iv) a(k) = 0 for all i < k < j. Based on numerical evidence obtained from 
high-precision computations for many different values of N (see Figure 4c for N = 161) we arrive at 

Conjecture 2: For 0 < n < IV — 1 the vector <fi n has exactly n zero crossings. 
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Figure 3: Comparing the first twelve eigenvectors <f> n (k) (black circles and red squares) with the scaled 
Hermite functions e 1 / 2, 0 n (efc) (black and red lines). Here N = 401, e = yj2n/N and the x-axis represents 
the index k E [—60,60]. 


Figure 4c also illustrates the following result 

l(4> n ) = 2L + 2\n/A]+l, (20) 

which follows easily from (14) and the above construction of <fi n . In particular, this result implies that 
about 25% of all elements 4> n (k ) of the eigenbasis are equal to zero. In view of Theorem 1 the following 
question naturally arises: are there other sets of eigenvectors {^ n }o<n<jV-i of the DFT such that tp n 
satisfies (20) and has exactly n zero crossings? In other words, is the eigenbasis (p n unique in a similar 
sense in which the vectors u n and v n are unique? 

3 Proofs 

For a nonzero vector a E we denote by R( a; z ) the rational function 

R(a;z) := ^ a(k)z k . (21) 
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(a) Hermite functions c 1 / 2 ip n (ek) 


(b) Eigenvectors <j) n (k) 


(c) The sign of cj) n (k) 


Figure 4: Intensity plots of Hermite functions e 1 ' 2, if n (ek), the eigenvectors 4> n (k) and the sign of the 
eigenvectors 4> n {k). Here N = 161, e = y/2i t/N, the x-axis represents the index n € {0,1,, 160} and 
the y-axis represents the index k G {—80, —79,... ,80}. Graph (c) shows the values where 4> n (k) = 0 
(black pixels), cf> n (k ) > 0 (grey pixels) and <f> n (k) < 0 (white pixels). 


We also denote ( := exp(— 2ir\/N). It is clear that b = Ta if and only if b(k) = N 1 / 2 i?(a; ( k ) for all 
k G In- The following simple observation will play the main role in the proof of Theorem 1: 

The function R( a; z ) has exactly /(a) — 1 zeros in C \ {0} (counting their multiplicity ). (22) 

To establish this fact we define j = min{fc G In ■ a(k) ^ 0} and note that z~^R( a; z) is a polynomial of 
degree 1(a) — 1. 

To give a preview of our method of the proof of Theorem 1, let us prove Kong’s result (8) (see [7, 
Theorem 1]) using the above property (22). Assume that for some nonzero vector a G C N we have 
1(a) — m and /(b) < N — m where b = Ta. Since /(b) < N — m it follows that b(k) = 0 for at least 
m indices k G In- In view of identity b(k) = A r_1 / 2 /?(a; ( k ), we see that the nonzero rational function 
R(a; z) has at least m distinct zeros on the unit circle, which is impossible because of (22). 

Proof of Theorem 1: Everywhere in this proof we will denote b = Ta. When we say that a vector or 
a function is “unique” we will mean “unique, up to multiplication by a nonzero constant”. 

Let us first prove part (i) for even vectors. Assume that N = 2 M + 1 for some M > 1, so that 
In = {— M, —M + 1,..., M}. It is clear that there exists a unique even vector a G R jV of length one 
(its only nonzero element is a(0)). The corresponding DFT vector b has elements b(k) = a(0)A^^ 1 / 2 
and has length N. Having dealt with this special case, we now assume that 1(a) = 2m + 1 for some 
1 < m < M — 1. If we impose the constraint 1(a) + /(b) = N + 1, then /(b) = 2 M — 2m + 1 and we 
must have R( a; ( k ) = b(k) = 0 for M — m + 1 < \k\ < M. This gives us 2m zeros of R( a; z), and since 
z m R( a; z) is a polynomial of degree 2m we can uniquely identify 

M 

R(^z) = z~ m (z-( k )(z-C l ), 

k=M—m+1 

which in turn uniquely identifies the vector a. Thus we have shown that for every 1 < m < M — 1 there 
exists a unique even vector a of length 2m + 1 satisfying 1(a) + /(b) — N + 1. Since there is also exactly 
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one even vector of length one and exactly one even vector of length N satisfying the same property (these 
were discussed above), we conclude that the total number of such vectors is M + 1. 

Consider next an odd vector a e M. N of length 2m + 1. Note that there are no odd vectors of length 
one, thus we can assume that 1 < m < M. Since a is an odd vector, the function f?(a; z) can be written 
in the form 

m 

R( a; z) = ^ a(k)(z k — z~ k ). 

k =1 

In particular, we see that R( a; ±1) = 0. Assuming that Z(a) + Z(b) = A^+3 we obtain Z(b) = 2M — 2m+3. 
Since b is an odd vector of length 2 M — 2m + 3 we must have R( a; ( k ) = b(k) = 0 for k = 0 and 
M — m + 2 < \k\ < M. This gives us 2m — 1 distinct zeros of R(sl;z). Taking into account the zero at 
z = —1 and using the fact that z m R(&;z ) is a polynomial of degree 2m, we can uniquely identify 

M 

R( a; Z) = z~ m (z - l)(z + 1) n (z - c k )(z - C~ k ). 

k=M—m +2 

Thus we have proved that for every 1 < m < M there exists a unique odd vector a of length 2m + 1 
satisfying /(a) + l(J- a) = N + 3. From our proof it should be clear that there can not exist an odd vector 
satisfying /(a) + l(J r a) — N + 1, because in this case the function R( a; z) would have 2m + 2 zeros, which 
is impossible due to (21). This completes the proof of part (i). 

Next, let us prove part (ii) for even vectors. Assume that N = 2M and that a e M. N is a nonzero 
even vector. In this case 


M—l 

R( a; z) = a(0) + ^ a(k){z k + z~ k ) + a(M)z M , 

k =1 

and assuming that a(M) = 0 we see that this function satisfies R( a; z ) = R( a; 1 /z). This identity implies 

= 0. (23) 

z=± 1 

There is a unique even vector of length one satisfying a(M ) = 0 (its only nonzero elements is a(0)), 
however for this vector we have b(M) = a(0)A r_1 / 2 ^ 0, thus we have to exclude this vector from 
consideration and we can assume that /(a) = 2m + 1 with 1 < m < M — 1. If /(a) + /(b) = + 2 then 

/(b) = 2 M — 2m + 1. Since b is also an even vector with b(M ) = 0, the equality /(b) = 2 M — 2m + 1 
implies that R{ a; ( k ) = b{k) = 0 for k = M and for M — m + 1 < \k\ < M — 1. Note that ( M = —1, 
thus R( a; z) has a zero at z — — 1, and from condition (23) we see that this zero must have multiplicity 
two. Thus we have determined all 2m zeros of the polynomial z m R( a; z ) of degree 2m, which allows us 
to identify uniquely 

M—l 

-R(a; z) = z~ m (z + l) 2 H ( z -( k )(z-(- k ). 

k=M—m+l 

To summarize: we have proved that for every 1 < m < M — 1 there exists a unique even vector a of 
length 2 m + 1 satisfying a(M) = b(M ) = 0 and /(a) + /(b) = A^ + 2. From our proof it is clear that 
there can not exist an even vector with the same properties satisfying /(a) + /(b) = N + 1, because in 
this case the function R( a; z) would have 2m + 1 zeros, which is impossible due to (22). This completes 
the proof of part (ii) for even vectors. 
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The proof of part (ii) for odd vectors follows the same steps and is left to the reader. 


□ 


The proof of Theorem 2 (and the derivation of the corresponding formulas in the Appendix A) is 
based on the following lemma. We emphasize that this result is true for all N > 1 (not only those 
satisfying N = l(mod 4)). 

Lemma 1. For all N >1 and m,n G {0,1,..., N — 1} we have 

N-l-n 

s ( n + k ) s ( N ~ k ~ l) e -^ (2fc+n+1 - 7V)(2m - n) = S(n)S(N + m-n- 1)S(N -m- 1). (24) 

fc =0 

Proof. Let us define the g-Pochhammer symbol as 

j - 1 

(a; q)j = na _ aqk )’ for J - 

fc =0 

and (z; q)o = 1. We begin with the g-Binomial Theorem 


(zq] q)n 


E 


(q;q) t 


„^ 0 (q; q)m(q; q)n- r 


-q 


m(m+l)/2 


l) m x m 


When x = q k , the above equality gives us 


(q;q)n+k 

(q; q)k 


n 


(q',q)r 




'.{Qi q)n—r 


g m ( m +l)/ 2 (^X) 


m km 

H ? 


(25) 


where we have also used the identity (g fc+1 ; q) n = (g; g) n+ fc/(g; g)*,. Next, we set g = exp(27ri/./V) in (25), 
we check that this choice of g implies 


(q-q) k = S{k)e- v ' lk{N - k ~ 1)/{2N \ for k > 0, 


and we obtain 


S(n + k) -mn(N-2k-n-l)/(2N) 

S(k) 


E 

m =0 


S(n) 


S(m)S(n — m) 


^ ^7rim(n-\-l) / N 


g2mkm/N 


The expression in the right-hand side is the inverse Fourier transform of the sequence given in square 
brackets. Applying the discrete Fourier transform results in 


S(n) 


-Qmg7rim(n-(-l)/ N yy—1 ^ ^ S(n T A.’) 2k—n—1)/(2N)^—2-irikm/N 


S(7n)S(n — m) 


k =0 


S(k) 


Formula (24) follows by simplifying the above expression and applying the identity (7). 


□ 
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Proof of Theorem 2: In order to prove formula (10), we set N = 4L + 1 and replace variables k, m 
and n in the equation (24) with fct —> L — n + k, m > L + n + l and n K > 2L + 2 n. This gives us the 
following identity 


L—n 

J2 S(3L + n + k)S(3L + n- k)e~ 2nikl/N = S(2L + 2n)S(3L -n + l)S(3L - n- l), 

k=—L-\-n 

which is valid for — L + n < k < L — n, —L — n < l < 3L — n and — L < n < L. Extending the above 
identity by periodicity for the values of / in the range —2 L < l <2L gives us the desired result (10). 

The fact that these vectors are the same ones that were discussed in Theorem l(i) follows from (10) 
and the fact that Z(u n ) = 2L — 2n + 1. In order to prove that these vectors are linearly independent, let 
us assume to the contrary that 


cquji + a 2 u i2 H-b a k u ik = 0, (26) 

where 0 denotes the zero vector, a 3 ^ 0 for all 1 < j < k and — L < i\ < i 2 < ■ ■ ■ < i k < L. Note that, 

by construction, u n {k ) = 0 for L — n + 1 < \k\ <2L and u n (k ) > 0 for \k\ < L — n. Therefore we have 

Ui k (L — if.) > 0 while Ui^L — i k ) = 0 for all other indices ii. Equation (26) then forces a k to be zero and 
we arrive at a contradiction. 

To prove formula (12), we use (9) and and check that 

4sin(7r(3L + n + k)/N ) sin(7r(3L + n + 1 + k)/N)u n -i(k) = u n (k + 1), 

4sin(7r(3L + n — k)/N ) sin(7r(3L + n + 1 — k) / N)u n -\{k) = u n (k — 1). 

Subtracting the above two equations and simplifying the result gives us formula (15). Applying the 
discrete Fourier transform to both sides of (15) we obtain 

2 L 

- 4sin(7r(2L + 2n)/N)(J r v n _ 1 )(l) = ^ e~ 2mkl/N (u n (k + 1) - u n (k - 1)) (27) 

k=—2L 
2 L 

= 2ism(2nl/N) e~ 2mkl ^ N u n (k) = 2isin(27rZ/A^)(J r u„)(Z). 

k=—2L 

In deriving the above identity we have used the fact that u n (±2L) = 0 for n > —L. Formula (12) follows 
from (10) and (27). The linear independence of the vectors v n can be established by the same argument 
as we used for u„. □ 


The following result will be needed for the proof of Theorem 3. We recall the definition of Catalan’s 
constant 


f7r/2 


G ■= 


ln(2sin(t/2))dt = 0.915965594... 


Lemma 2. There exists an absolute constant C > 0 such that for all N large enough and all k G 
[- N 2 /3, N 2 / 3 ] n z we have 

u 0 (k) = 2- 1 / i Ne NG/n - nk2/N {l + S) (28) 

and \£\ < CN- 1 / 3 . 
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Proof. First of all, we use property (7) and check that 


u 0 (k) = S(3L + k)S(3L - k) = N (29) 

Let us denote f(x ) := ln(2 sin(nx/N)) for 0 < x < N. Applying the Eulcr-Maclaurin summation formula 
we obtain 


ln(S(3L + k)/S(L + k)) 


?> L + k j>3L+k 

J2 M= f(x)dx+(f(3L + k)~f(L + k))/2 

T. 4-t-l-l J L+k 


j=L+k+l 


+ (f(3L + k)-f(L + k))/12 + H, 


(30) 


where the remainder term has an upper bound 


W <55 


f>3 L+k 


\f"(x)\dx. 


'L+k 


We compute the first term in the right-hand side in (30) 


r»3 L+k 


h ■= I fix) dx = -^(C1 2 (2 tt(3L + k)/N) - Cl 2 (2vr [L + k)/N)), 


I L+k 


where the Clausen function is defined by 


Cl 2 (t) := — / In |2sin(u/2)|dw. 


The function Cl 2 (f) is odd and periodic with period 27T (see [8, Chapter 4]), which gives us 

N 

h = —(C1 2 (2t t(L + 1 - k)/N) + C1 2 (2t t(L + k)/N)). 

Zn 


The Clausen function can be written as Taylor series C1 2 (7 t/2 + z) — ^A >0 a 3 zf where ao = G, a± = 
— ln(2)/2, a 2 = —1/4 and the series converges if \z\ < 7t/ 2. Using this result we obtain 

h = J2 a ji 27l /N) j ~\i3/4: - k) j + (k- l/4) j ) = NG/n - ln(2)/4 - -^(2 k 2 -2 k + 5/8) 
j> o 

+ (A)' (3fc 2 /2 - 3k/2 + 13/32) + £ % (2 I r/JV)^ 1 ((3/4 - k)> + (k - 1/4)+ 

j> 4 


Assuming that \k\ < N 2 / 3 we have 


h = NG/tt - ln(2)/4 - nk 2 /N + 0(N~ 1/3 ) + Oj(2 tt/JV) J ' _1 ((3/4 - k) j + {k - 1/4 )>). 

j> 4 

When N > (47t) 3 , each term in the series in the right-hand side of the above formula can be estimated 
as follows 


|(27r/A0 J_1 ((3/4 - k) j + (k - 1/4) J ')| < 2(N 2 / 3 + l) j (2n) j ~ 1 N 1 - j 

< (4tt yN 1 -^ 3 < {a^Yn- 1 / 3 , 
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for all j > 4. Combining all the above two formulas and using the fact that the series Ylj>o \ a i\ converges 
we arrive at our final estimate of Ji, which is valid for all N > (47t) 3 and all \k\ < N 2//3 : 

h = NG/n - ln(2)/4 - vr k 2 /N + 0(iV 1/3 ). 


Dealing with the remaining terms in the right-hand side of (30) is much easier. We check that the 
second term satisfies 


I 2 ■.= (f(3L + k) - f(L + k))/2 = - In 


( sin(^ + 7 t(3/4 - k)/N)\ 


Vsin(| + n{k - 1/4)/IV) J 
Similarly, the third term can be estimated as 


= - ln(l + 0(N~ 1/3 )) = 0(N~ 1/3 ) 


h '■= ( f'(3L + k) - f'(L + k ))/12 = p^(cot(7r/4 + 7 t(3/4 - k)/N) + cot(7r/4 + n(k - l/4)/lV)) = 0(N 3 ) 
Finally, the remainder can be bounded as follows 


\n\ < — 

1 1 12 


f*SL-\-k 


7T 


r»3 L+k 


/ \f"(x)\dx — 

I L+ k W 12N 2 


csc(7Tx/N) 2 dx 


I L+k 


< 


7l~ 


12N 2 


x (N/2) x max {esc (7ra; / A^) 2 : N/ 4 - N 2/3 <x< 31V/4 + N 2/3 } = OiN- 1 ). 


Combining (29) with (30) and the above four asymptotic results ends the proof of (28). 


□ 


Proof of Theorem 3: We recall that e = s/2n~JN and we define rj := 2 1 / 4 A^exp(A^G/7r) and 

J N :={keZ : \ek\ < 21V 1 / 16 }. 


Our first goal is to prove the following upper bound: there exists an absolute constant C\ > 0 such 
that 

0 < u n (k) < Ci? 7 exp(—IV 1 / 8 ), for \n\ < 2 and k E In\ Jn■ (31) 

In order to establish (31) we note that 


u n (k + 1) 


sin(7r(L — n — k)/N) 

— - — ' j : T+^ u n{k), 

sm(7r(L — n + k + 1)) 


which shows that u n (k + 1) < u n {k ) for k > 0. By symmetry we also have u n (k — 1) < u n (k ) for k < 0. 
Lemma 2 implies that when ek is close to 21V 1 / 16 then Uo(k) < 2^exp(—21V 1 / 8 ). Since u n (k ) decreases 
as \k\ increases, this bound holds true for all k G In \ Jn, and we have proved (31) for n — 0. For other 
values of n we use the same method coupled with the identity 


u n+ i(k) — 2 (cos(2nk/N) — cos(2n(L — n)/N))u n (k), (32) 

which follows at once from (13). 

Next, we introduce t k := 2 sin(e 2 /c/2)/e. Note that as e —» 0 we have t k = ek + 0(e 5 k 3 ). Using the 
double-angle formula for the cosine function we rewrite (32) in the form 

u n +i(k) = (2 - 2sin(e 2 (?r + 1/4)) - e 2 t 2 k )u n {k) 

= (2 - e 2 (t 2 + 2n + 1/2) + e 6 (n + l/4) 3 /3 + 0(e w ))u n (k). (33) 
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We emphasize that the implied constant in the error term 0(e 10 ) does not depend on k. Applying (33) 
we obtain 


u 0 (k) = (2 - e 2 (t 2 - 3/2) + 0(e 6 ))u-i(k), (34) 

«i(*0 = (4 - e 2 (44 - 2) + e 4 (4 - t 2 k - 3/4) + 0(e 6 ))u^(k). 

The identity S(2L) = \/iV coupled with (5) give us 

N~ 1/2 S(2L + 1) = 2sin(7r(2L + 1)/N) = 2cos(e 2 /4) = 2 + 0(e 4 ), (35) 

N~ 1/2 S(2L + 2) = 4sin(7r(2L + 1)/N) sin(7r(2L + 2)/iV) = 4cos(e 2 /4) cos(3e 2 /4) = 4 - 5e 4 /4 + 0(e 8 ). 

Using formulas (34) and (35) and the fact that v n {k) = sm(2nk / N)u n (k) we arrive at the following result 

w 0 (k) = 2 u 0 (k) = (4 - e 2 (2 t 2 k - 3) + 0(e 6 ))n_i(/c), 

x 0 (k) = v 0 (k ) + N~ l t 2 S(2L + 1)d_ 1 (/c) = (4 + e 2 (3/2 — i£) + 0(e 4 )) sin(27rfc/AT)it_i(A;), 

3/o(fc) = -ui(/;) + N~ 1/2 S(2L + 2)it_i(fc) = (e 2 (4^ - 2) - e A (t\ -t 2 k + 1/2) + 0(e 6 ))n_i(fc), (36) 

Zo(k) = —v 0 (k) + N~ 1 ^ 2 S(2L + l)n_i(fc) = (e 2 (t 2 , — 3/2) + 0(e 4 )) shi(2nk/N)u_i(k). 

Again, in all of the above formulas the implied constants in the error terms 0(e J ) do not depend on k. 
Now we restrict k to the interval Jn■ For k e Jn we have \t k \ < AN 1 / 16 , thus et k = 0(N~ l t A ). We also 
check that sm(2nk/N) = etk + 0((e4) 3 )) and 

u-i{k) = | exp(—(efc) 2 /2)(l + 0(e 2/3 )), k e J N , (37) 

which follows from (33) and Lemma 2. Combining all of the above results with the following formulas 
for Hermite polynomials 

H 0 (x) = 1, H x (x) — 2x , H 2 (x ) = Ax 2 — 2, H 3 (x) = 8x 3 — I2x, 

we conclude that for k G Jn 

w 0 (k) = 2r]H 0 (t k ) exp(—(efc) 2 /2)(l + 0(e 2/3 )), 
x 0 (k) = erjH^tk) exp(-(e/c) 2 /2)(l + 0(e 2/3 )), 

Vo(k) = e 2 {j]/2)H 2 {t k ) exp(—(efc) 2 /2)(l + 0(e 2/3 )), (38) 

^o(fc) = e 3 (r)/16)H 3 (t k ) exp(-(efc) 2 /2)(l + 0(e 2/3 )). 

Let us consider, for example, vector z 0 . We use the above result and and the estimate (31) to conclude 

256e _5 7U 2 ||z 0 || 2 = 256e _5 ?U 2 ^ £ 0 (&) 2 + e ^ i7 3 (4) exp(-(efc) 2 /2)(l + 0(e 2/3 )) 

k£lN\JN AiGJat 

= 0(AT 7 / 2 exp(-2AT 1 / 8 )) + e ^ if 3 (t fc ) exp(-(efc) 2 /2)(l + 0(e 2/3 )). 

&G JW 

Thus, as iV —>■ +cx) we have 
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Recall that 0 3 = zo/||zo||. Then (39) combined with the last formula in (38) show that for all k G Jn 
we have 

e" 1/2 0 3 (/c) = ^ 3 (efc)( 1 + o(l)). 

In the range k G In\Jn both the left hand side and the right-hand side are bounded by 0(N A exp(—IV 1 / 8 )) 
(see (31)), thus as N —> +oo both of these expressions converge to zero, uniformly in k G In \ Jn- This 
concludes the proof of the convergence result (18) for n = 3. The proof in the cases n G {0,1, 2} follows 
exactly the same steps. 

The proof for the remaining cases n G {4, 5,6,7} will require the following simple fact, whose proof 
is left to the reader. 

Fact: Assume that the orthonormal vectors c and d are obtained from two linearly independent vectors 
a and b by Gramm-Schmidt orthogonalization 


c = 



d 



d 



(40) 


Then the vectors c and d would be the same if instead of a and b we started with a = ^a and b = 
/3 2 b + /3 3 a for some ( 3 1 > 0 and /3 2 > 0. 

Now, armed with this result, let us prove the convergence result (18) for n = 4. According to the 
Definition 1, the vectors <f> 0 and 0 4 are obtained by the Gramm-Schmidt orthogonalization (40) starting 
from vectors w 0 and wq . According to the above fact, we will obtain the same result if we start with 
vectors 

w 0 := (2^) _1 w 0 , and wq := (32/r])e~ 4 (wq — (2 — e 2 + e 4 /16)w 0 ) . (41) 

It is clear from (38) that 

w 0 (k) = H 0 (t k ) exp(—(e/c) 2 /2)(l + 0(e 2/3 )) (42) 

for k G Jn- Using (34) and (35) we compute 


Wi(k) = Ui(k) + N 1 ^ 2 S(2L + 2)u-i(k) = (8 — e 2 (4f(fc) 2 — 2) + e A (t(k) 2 — t(k ) 2 — 2) + 0(e 6 ))w_i(A:). 


Combining this result with (36), (37) and the following expression for the Hermite polynomial 

H±(x) = 16x 4 — 48x 2 + 12, 


we find 


wi(k) = (2 /rj) x (16 1\ - 48 t 2 k + 12 + £)u_i(k) = (. H 4 (t k ) + S) exp(-(efc) 2 /2)(l + 0(e 2/3 )), (43) 

where the error term satisfies £ = 0(e 2 (l + f 4 )). For k G Jn we can estimate £ = O(e). Thus we have 
proved that for all k G Jn 

wi{k) = H 4 (t k ) exp(—(efc) 2 /2)(l + 0(e 2/3 )). (44) 

It is clear from (31) and (41) that for k E In\ Jn we have an upper bound l'uq(fc)! = 0(N 2 exp(— N 1 ^ 8 ), 
thus w\ (k) converges to zero as N —>■ +oo uniformly in k G In \ Jn- 
Now, the vector <f> 4 is given by = e/||e||, where we set 


W| • Wn A 
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Applying (42), (44) we conclude (in the same way as we did above for ||z 0 || 2 ) that as N —> +oo 


e w o = e 


: ^ w 0 (k) 2 —> / exp(— x 2 )dx = \Jn, 

kei N ,, 


ewi ■ w 0 = e wo(k)wi(k ) —> J H 4 (x) exp(— x 2 )dx = 0. 

fceijv j* 


Thus the quantity wy • w 0 /||w 0 || 2 converges to zero as N —)■ +oo, which shows that e(k) —> Wi(k) as 
N —>- Too (uniformly in k G In)- This result, combined with (44) proves that the convergence result 
(18) holds true for n = 4. 

The proof of (18) for n G {5, 6, 7} follows the same steps, though the computations become more te¬ 
dious (the use of a symbolic computation package is highly recommended for checking these expressions). 

For example, when n — 7 we use (16) and (33) to find 

z 0 (k ) = (e 2 (2 t 2 k — 3) — e 4 (4 — 5 t 2 k T 43/8) T e 6 (t 2 + 1)/16 T 0(e 8 )) sm(27rk/N)u- 2 (k), 

Zi(k ) = (e 2 (12f 2 , — 18) — e 4 (6f^ — 18 t 2 k + 57/4) T e 6 (f® — 9f^/2 T 11^/4 T 165/4) T 0(e 8 )) sin(2nk/N)u_ 2 (k). 

We define zi := (512/r7)e^ 7 (zi — (6 — 6e 2 T 49e 4 /16)z 0 ) and check that 

z\{k) = (4/r7)e _1 (128t®, — 1344t/. T 3360 1 2 — 1680 T £) sm(2irk/N)u- 2 (k) 

= ( H 7 (t k ) T £) exp(—(efc) 2 /2) cos(7rfc/A^)(l T 0(e 2/3 )), 

where \£\ < e 2 (l T t\). Starting with this result and following exactly the same argument as we have 
used in the case n — 4 one can show that the convergence result (18) holds true for n = 7. The details 
and the proof of the remaining cases n G {5,6} are left to the reader. 
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Appendix A 

In this section we present the analogues of Theorem 2 and Corollary 1 for the remaining cases N G 
{0,2,3} (mod 4). In each case we give formulas for the even and odd vectors u n and v n (which were 
described in Theorem 1) and compute their corresponding discrete Fourier transforms. The proof of these 
results is based on Lemma 1 and it is very similar to the proof of Theorem 2, we omit all the details. 
We construct a basis for each eigenspace W, X , Y and Z, from which one can build the Hermite-type 
eigenvectors cj) n in exactly the same way as in Definition 1. 

We recall that the interval In and the discrete Fourier transform T are defined by (2) and the 
sequence S(k) is defined by (5). 


A.l The case when N = 4L 

We define 2L+1 vectors {u„}_i<„<i and 2L—1 vectors {v n }_£ + i< n <x,_i as follows: for —2L+1 < k <2L 
we set U-L(k) := 2L and 

u n (k) := cos(nk/N)S(3L + n — 1 + k)S(3L + n — 1 — k), —L + 1 < n < L, 
v n (k) := sm(7rk/N)S(3L + n — 1 + k)S(3L + n — 1 — k), —L + 1 < n < L — 1. 


These N vectors are linearly independent and they satisfy 

, /9 S(2L + 2n) 

J-U n = N I ;:—7— — - TTTTV U _n, ~L < U < L — 1, 

-v_ n , — L -\-1 < n < L — 1, 


2 sin(7r(3L + n) /N) 
= —L/V -1 / 2 - g(2L + 2n) 


2 sin(7r(3L — n)/N ) 
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and Tu l = 2a//Vu_^. The vectors {u n }_£ + i< n <L + i and {v n }_L + i< n <£_i are precisely the vectors 
described in Theorem 1 (ii). An equivalent expression for the vectors u n is given by 


2L—1 

u n (k) = N2 L+n cos(nk/N) 2 (cos(27rfc/A) — cos(2nj /N)). 

j=L—n +1 

Define A vectors {w n } 0 <n<L, {x n } 0 < n <L-i, {y n }o<n<L-i, {z n }o<n<L -2 as follows: 

-i /9 S(2L + 2n) 

w n = u n + N A ■ f for , w A n u -" ’ 0 < n < L — 1, 
2sm(7r(3L + n)/N) 

w L = u L + 2A 1/2 u_ L , 

, /2 S'(2L + 2n) 

x n = v„ + A 1 -rxy v -n , 0 < n < L - 1, 

2 sm(7r(3L — n)/N) 

Ar , /9 SY2L + 2n + 2) 

y n = -u n+1 + A A . , 7 xy u _ n ~i , 0 < n < L — 2, 

2sm(7r(3L + n + lJ/A) 


y L -i = -u L 


+ 2A 1 / 2 u_ L , 


at _ 1/2 S(2L + 2n + 2) 

Zn — ~v n+ i + A A . . . -X7wv_ n _ 1 , 0 < n < L — 2. 

2 sm(7r(3L — n — 1) / A ) 

The vectors {w n }o<n<L (respectively, {x„}o<n<L-i, {yn}o<n<L-i and {z n }o<n<L- 2 ) form a basis for the 
eigenspace W (respectively, X, Y and Z). 


A.2 The case when N = 4L + 2 

We define 2L + 2 vectors {u n } -l<u<l+i and 2 L vectors {v n }-L+i<n<L as follows: for —2 L < k < 2L + 1 
we set u-L(k) := 2L + 1 and 

u n (k) := cos(nk/N)S( y 3L + n + k)S(3L + n — k ), —L + 1 < n < L + 1, 

n n (/c) := sin(7r/c/A)S'(3L + n + k)S(3L + n — k), —L + 1 < n < L. 

These A vectors are linearly independent and they satisfy 


Xu 


n 


A - 1 / 2 


S{2L + 2n) 

2sin(7r(3L + n + 1)/A) 1 n ’ 


—L < n < L, 


Tv 


n 


-iA " 1 / 2 


S{2L + 2n) 

2 sin(7r(3L — n + 2)/A) 1 


—L + 1 < n < L, 


and Ju L+ i = 2A 1 / 2 u_£. The vectors {u n }_£ + i< n <£ and {v n }-L+i< n <L are precisely the vectors de¬ 
scribed in Theorem 1 (ii). An equivalent expression for the vectors u n is given by 


2 L 

u n {k ) = N2 L+n cos(nk/N) 2 (cos(27t/c/A) — cos(27tj’/A)) . 

j=L—n +2 
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Define N vectors {w„}o<n<L, {x n }o<n<L-l, {yn}(Kn<L, {z n }o<n<L-l as follows: 


, /2 S(2L + 2n + 2) 

Wn = u n+ i + N / . 7 ( u_ n , 0 < n < L - 1, 

2 sin(7r(3L + n + 2) /A) 


w L = U L+ i 


+ 2iY- 1/2 u_ L , 


at-i /2 S(2L + 2n + 2) 

X n = V n+1 + N ' --ww V_ n , 0 < n < L - 1, 

2 sm(7r(3L — n + 1)/iv) 

, /2 S'(2L + 2?i + 2) 

yn = -Un+1 + iV 1 ———— , 9 w An U-n , 0 < 71 < L - 1, 

2 sm(7r(3L + 71 + 2)/N) 


y l = -u L+1 


+ 2iV- 1 / 2 u_ i , 


, r _ 1/2 S(2L + 2n + 2) 

z n = -v n+1 + N 1 - lW , n v -n , o < n < L — 1. 

2 sin( 7 r( 3 L — n + l)/iv) 


The vectors {w n } 0 < n <L (respectively, {x n } 0 < n <L_i, {y„}o<n<L and {z ra } 0 < n <L-i) form a basis for the 
eigenspace W (respectively, X, Y and Z). 


A.3 The case when N = 4L + 3 

We define 2L+2 vectors {u„}_z / < n </ /+ i and 2L+1 vectors {v n }-L<n<L as follows: for —2L—1 < k < 2L+1 
we set 

u n (k) := S(3L + n + 1 + k)S(3L + n + 1 — k), —L < n < L + 1, 
v n (k) := sm(2iTk/N)S(3L + n + 1 + k)S{3L + n + 1 — k), — L < n < L. 

These N vectors are linearly independent and they satisfy 

T\i n = N~ 1 ^ 2 S(2L + 2n)ui_ n , — L < n < L + 1 
Jv„ = — S{2L + 27i + l)v_„, — L < n < L. 

The vectors u„ and v n are precisely the vectors described in Theorem l(i). An equivalent expression for 
the vectors u„ is given by 


2 L 

u n (k) = N2 L+n (cos(2nk/N) — cos(2nj/N)). (45) 

j=L—n +2 

Define N vectors {w n } 0 < n <L, {x n } 0 <n<L, {y n }o<n<L, {z n }o<n<L-i as follows: 

w n = u n+ i + N~ 1/2 S(2L + 2n + 2)u_ n , 0 < n < L, 
x„ = v n + N~ 1/2 S(2L + 2 n + l)v_ n , 0 < n < L, 
y n = -u n+1 + N~ l/2 S(2L + 2n + 2)u_ n , 0 < n < L, 
z n = —v n+1 + N~ 1 / 2 S(2L + 2n + 3)v_ n _! , 0 < n < L — 1. 

The vectors {w n } 0 < ri <z / (respectively, {x n } 0 < n < L , {y„} 0 <n<L and {z n } 0 < n <i_i) form a basis for the 
eigenspace W (respectively, A", Y and Z). 
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